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[1]  INTRODUCTION 

In  order  to  accomplish  work  in  any  rotor,  angular  momentum  must  be 

Imparted  to  the  fluid.  The  lifting  action  of  the  rotor  blades  produces 

cross  channel  pressure  gradients.  At  a sufficient  distance  from  the  inner 

wall,  viscous  effects  are  negligible  and  the  pressure  gradients  are 

balanced  by  streamline  curvature.  Close  to  the  inner  wall,  a boundary  layer 

exists.  The  tiuid  within  this  boundary  layer  does  not  have  sufficient 

momentum  to  balance  the  pressure  gradients  imposed  by  the  inviscid  outer 

flow.  The  result  is  a cross-flow  component  containing  vorticity  aligned 

in  the  streamwise  direction.  This  additional  streamwise  vorticity 

causes  a deviation  in  the  rotor  outlet  angles.  In  some  cases,  this 

streamwise  vorticity  causes  a roll-up  of  the  low  momentum  fluid  near 

* 

the  wall  into  a vortex-type  flow  [1]  . 

There  are  no  simple  techniques  that  accurately  predict  these  inner 
wall  velocity  gradient  effects  for  rotational  flows.  It  is  important 
to  recognize  that  a satisfactory  description  of  this  boundary  layer  flow 
cannot  come  from  refinements  of  the  two-dimensional  boundary  layer  theory, 
but  rather  from  three-dimensional  analysis.  The  reason  is  that  boundary 
layer  behavior  in  rotors  exhibit  variations  in  its  lateral  direction.  The 
occurrence  of  these  secondary  flows  which  are  boundary  layer  flows  having 
a component  normal  to  the  mainstream  direction  arise  principally  from 
blade  to  blade  and  radial  pressure  gradients.  Similarity  between  the 
results  of  secondary  flow  analysis  and  the  cross  flow  in  the  outer  part 
of  a three-dimensional  boundary  layer  has  been  established  by  Horlock  [2]. 

The  distributed  passage  secondary  vorticity  generated  by  turning  a boundary 

layer  flow  is  one  of  three  possible  types  of  secondary  streamwise  vorticity, 

_ 

Numbers  in  brackets  refer  to  documents  in  the  list  of  references. 
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which  can  occur  near  the  inner  of  a rotor.  As  discussed  by  Hawthorne  [3], 
three  components  of  secondary  vorticity  are  identified  in  the  direction 
of  flow  at  the  exit  of  a blade  row.  These  secondary  sources  of  vorticity 
are  usually  regarded  as  a perturbation  on  the  primary  flow.  They  are 
(1)  a distributed  passage  vorticity  in  the  blade  passage  which  may  result 
in  the  formation  of  the  so-called  passage  vortex,  (2)  the  trailing  shed 
vorticity,  and  (3)  the  trailing  filament  vorticity.  The  latter  two  types 
of  vorticity  are  due  to  the  vortex  sheet  leaving  the  blade  trailing  edge 
and  lead  to  the  formation  of  another  vortex  which  is  opposite  in  rotation 
to  that  of  the  passage  vortex  as  shown  in  Figure  1.  The  trailing  filament 
vorticity  is  caused  by  the  stretching  of  the  vortex  filaments  as  they 
move  over  the  surfaces  of  the  blades.  The  last  component  of  vorticity 
is  the  trailing  shed  vorticity  which  is  caused  by  the  variation  of  circulation 
along  the  span  of  the  blades. 

With  the  exception  of  trailing  shed  circulation  which  exists  along 
the  blade  even  when  the  incoming  flow  is  uniform,  each  of  these  vorticity 
components  could  be  attributed  to  the  existence  of  the  wall  boundary 
layer  and  not  the  change  of  the  boundary  layer  due  to  viscous  effects 
as  flow  passes  through  the  rotor.  Therefore,  the  primary  assumption 
leading  to  the  existing  theoretical  descriptions  of  secondary  flows  is 


that  viscous  effects  produce  a boundary  layer  on  the  wall  upstream  of 
the  rotor.  Whereas  within  the  rotor,  the  imposed  pressure  gradients 
play  the  major  role  and  viscosity  has  little  effect  on  the  resulting 
secondary  flows.  This  assumption  is  characteristic  of  what  is  generally 
termed  inviscid  secondary  flow  analysis. 

As  shown  in  an  analysis  of  secondary  flows  by  Came  and  Marsh  [4], 
the  total  strength  of  these  three  secondary  components  of  streamwise 
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vorticlty  is  zero  in  the  flow  downstream  of  a many  bladed  cascade. 

However,  this  vorticity  in  the  flow  does  have  an  effect  on  the  flow  field. 
The  primary  effect  of  this  secondary  vorticity  is  a deviation  in  the 
blade  outlet  angle  due  to  the  passage  secondary  vorticity.  This  deviation 
can  be  quite  large  when  the  incoming  velocity  gradient  to  the  rotor  is 
large  or  when  the  flow  is  turned  through  a large  angle.  The  solution  for 
this  deviation  must  be  consistent  with  the  trailing  vortex  sheet,  but  the 
strength  of  this  sheet  does  not  need  to  be  known. 

Another  effect  of  secondary  distributed  vorticity  generated  near 
the  wall  of  a rotor  is  evidenced  by  the  existence  of  a cavitating  vortex. 
The  appearance  of  the  cavitating  vortex  is  similar  to  that  shown  in 
Figure  2.  The  structure  of  the  vortex  varies  with  velocity  gradient-rotor 
configuration.  The  magnitude  of  the  minimum  pressure  coefficient  (Cpmin) 
associated  with  this  vortex  depends  on  both  the  primary  and  secondary 
vorticity.  However,  the  resultant  strongly  swirling  flow  appears  to  be 
organized  by  the  secondary  vorticity. 

Theoretical  understanding  of  secondary  flows  through  a rotor  is 
aided  by  approximate  solutions  of  the  fluid  flow  equations  that  govern 
the  flow  process.  The  equations  themselves  are  approximate  because 
certain  simplifying  assumptions  must  be  made  before  the  solution  can  be 
obtained.  However,  any  secondary  flow  theory  which  is  applicable  to  the 
vortex  problem  must  be  least  include  the  effects  of  blade  twist,  variable 
relative  velocity  through  the  rotor,  large  velocity  gradients  and 
rotation. 


The  objective  of  this  report  is  to  review  and  discuss  secondary  flow 
theories  which  can  be  employed  to  estimate  the  vorticity  created  near  the 
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inner  wall  of  a rotor.  There  are  many  different  approaches  to  the  derivation 
of  secondary  flow  equations;  however,  Horlock  [5]  in  a discussion  of  a 
paper  on  secondary  vorticity  in  axial  compressor  blade  rows  shows  that 
the  different  approaches  all  lead  to  essentially  the  same  result.  A 
helpful  review  of  secondary  flow  theories  can  be  found  in  the  papers  of 
Lakshminarayana  and  Horlock  [6],  the  staff  of  NASA  Lewis  Research  Center 
[7],  Hawthorne  and  Novak  [8],  Lakshminarayana  and  Horlock  [9]  and  Salvage 
[10]. 

[2]  PASSAGE  SECONDARY  VORTICITY  IN  A ROTOR 

The  equation  of  motion  for  incompressible  flow  with  reference  to  axis 
rotating  at  constant  angular  velocity  (0)  is  given  in  Greenspan  [11]  as 

fxW  + itixV  = - V’(^)  - V'[W2/2  - l/2(fixr)  • (ftx£)]  - (u/p)V,xV'xW  (1) 

where  the  prime  denotes  differentiation  with  respect  to  the  rotating  frame 
(fl) , W is  the  relative  velocity  and  £ is  the  relative  vorticity  given  by, 
£“V'xW=u)-2 ft.  Taking  the  curl  of  Equation  (1)  gives  the  vorticity  equation 

(W-V')f  = (f*V')W  - 2V'x(fixW)  + (y/p)V,2£  . (2) 


Expessions  for  the  absolute  secondary  vorticity  (w)  defined  along  a 
relative  streamline  valid  for  a rotor  are  obtained  from  Equation  (2)  by 
Lakshminarayana  and  Horlock  [9].  The  resultant  equations  are 


2lV  2nb,(V 

WR'  w2 


— + (y/p)s? 


(3) 
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3n' 


[wn,W] 


WU)^ 

T 


WV  3ab' 

ab,  3s' 


+ (y/p)n 


V'2? 


(4) 


where  as  shown  in  Figure  3,  s',  n' , b'  represent  the  natural  coordinates  for 
the  relative  flow,  W is  the  relative  velocity,  and  are  absolute 

vorticity  resolved  along  the  relative  streamline  (s')  and  the  principal 
normal  direction  (n'),  ft  is  the  rotation  vector  and  R'  is  the  radius  of 
curvature  of  the  relative  streamline.  In  this  relative  coordinate  system 
indicated  by  the  primes,  the  streamwise  direction  (s')  is  defined  along  the 
flow  direction,  the  principal  normal  direction  (n')  is  defined  toward  the 
center  of  curvature  of  the  steamlines,  and  the  bi-normal  direction  (b') 
is  defined  as  s'x  n’=b'.  Further  definitions  and  descriptions  of  the  flow 
in  this  relative  coordinate  system  are  given  by  Lakshminarayana  and  Horlock  [9]. 

The  means  by  which  the  streamwise  component  of  vorticity  is  produced  in 
this  relative  flow  are  similar  to  those  in  a stationary  system.  However,  it 
is  important  to  note  that  additional  secondary  vorticity  is  generated  when 
ftxW  has  a component  in  the  relative  streamwise  direction.  Rotation  has  no 
effect  when  the  absolute  vorticity  vector  lies  in  the  s'-n'  plane  and  the 
rotation  (ft)  has  no  component  in  the  bi-normal  direction  (b'). 

The  equations  for  secondary  vorticity  created  in  the  passage  namely 
Equations  (3)  and  (4)  are  extremely  difficult  to  apply  because  they  are 
partial  differential  equations.  Therefore,  some  simplifying  assumptions 
and  specification  of  the  flow  field  must  be  made  for  each  application.  The 
basic  flow  parameters  for  this  case  discussed  in  this  paper  are  a thick 
boundary  layer  entering  the  rotor  and  the  rotor  operating  at  a flow 
coefficient  near  1.0.  It  is  also  important  to  note  that  the  primary 
flow  through  the  rotor  is  rotational. 


-10- 


17  August  1977 
MLB: jep 


» 


In  order  to  evaluate  the  relative  importance  of  the  terms  that  describe 
the  passage  secondary  vorticity,  the  vorticity  equations  were  first 
nondimensionalized.  Then  an  order  of  magnitude  evaluation  can  be  used  in 
the  selection  of  the  terms. 

The  normal  vorticity  equation  can  now  be  evaluated  in  the  following 
manner: 


V - 

w /6  V 

o 


V 


W /R' 
o m 


b f 


where 


W 


o 


6 


£ 

*R 


r 


o 


'v  initial  relative  velocity 
^ boundary  layer  thickness 
^ camber  radius  at  mid-span 

'V'  distance  between  neighboring  streamlines  in  b' 
direction 

blade  span 

^ rotor  tip  radius 

^ rotor  inner  radius. 


Applying  these  nondimensional  variables  to  Equation  (4)  and  ignoring 
viscous  terms  in  the  normal  direction,  we  find 
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or 
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w) 
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2 
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W 

0 

"V 
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o 

aIb' 

R’*r 

m 

T 

» 4 

6*c 

V 

(5-a) 


3 

37' 


(a)  , W} 

n 


6»c 

R'  -r 
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W (iJn,  3ab, 
- ■ • - ■ - - 

a,  , 3s' 

D 


(5-b) 


The  ratio  cf  boundary  layer  thickness  to  rotor  tip  radius  is  of  order  one 
for  the  applications  considered  in  this  paper.  Also,  the  blade  chord  (c) 
is  an  order  to  magnitude  less  than  the  camber  radius.  Therefore,  the 
normal  relative  vorticity  equation  becomes 


W)  = - 


W to  , 3a,  , 
n b 

V as' 


(6) 


Equation  (6)  can  be  integrated  by  combining  the  two  terms  into  one 
differential.  The  result  is 


W. 


0) 


Wn| 


1 3bl 

W V 


(7) 


where  the  subscript  1 refers  to  the  rotor  inlet  and  unnumbered  subscripts 

refer  to  any  position  along  a streamline  in  the  rotor.  This  simplified 

equation  shows  a dependence  of  the  relative  component  of  the  normal 

vorticity  on  the  variation  of  relative  velocity  (W)  along  a streamline, 

on  the  relative  component  of  the  absolute  inlet  normal  vorticity  (w^) 

and  on  the  flow  convergence-divergence  in  the  bi-normal  direction  (a  ,). 

Equation  (7)  can  be  evaluated  along  a relative  streamline  if  the 

streamline  spacing  (a,  .)  in  the  bi-normal  direction  is  known.  This  can 

b 

be  easily  estimated  by  applying  a numerical  technique,  such  as  the 
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streamline  curvature  method,  to  the  primary  flow  field.  For  some  rotors 
having  a high  hub  to  tip  ratio,  the  streamline  spacing  in  the  radial  direction 
(a^,)  remains  constant  and  Equation  (7)  then  becomes 


U) 


U)  • 

nl 


W 

_2 

w 


(8) 


This  result  relates  the  change  in  the  relative  component  of  the  absolute 
normal  vorticity  (u^,)  to  changes  in  the  relative  velocity  (W) . Also, 
assuming  a weak  shear  flow  so  that  the  axial  velocity  remains  constant 
gives 


to 


cos  8 
cos  8 


(9) 


The  last  two  relationships  are  commonly  used  for  secondary  flow  cascade 
theory. 

In  a similar  manner,  the  streamwise  secondary  vorticity  equation 
can  be  nondimensionalized  and  evaluated  as  follows: 


to 


c 

R' 


W 


to 


n. 


to  sin 


V 


H sin  $ 


n 


n 
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where 

4*  'V  rake  angle  of  the  rotor 

ft  'V  magnitude  of  rotational  vector  of  rotor. 

Applying  these  additional  nondimensional  variables  to  Fetation  (3) 
for  the  streamwise  vorticity  we  get 


6 , 

■ 

V 

2d)  , Rf  « sin  V 

n m 

► « 

2 V V 

’ * 

2 5n'  V 

3s’ 

W 

. d 

“ 1 w 

W R'  o 

w2 

. 4 

W2 

. 4 

4 — . 2_ 

6 Re 


-2  - 

1 3 “s’ 


W2  3 b’2 


(10) 


The  second  term  on  the  right-hand  side  of  Equation  (10)  is  of  order  one 
because  the  flow  coefficient  is  of  order  one.  This  term  is  due  to 
curvature  induced  secondary  vorticity.  The  third  term  is  also  of  order 
one  and  is  due  to  rotation- induced  secondary  vorticity.  However,  both 
of  these  terms  are  as  important  to  secondary  flow  development  as  the 
first  term. 

The  fourth  term  on  the  right-hand-side  of  Equation  (10)  is  a viscous 
dissipation  term  which  can  be  neglected  in  most  cases.  Even  if  the 
Reynolds  number  is  based  on  local  parameters  in  the  boundary  layer,  this 
term  is  still  small  compared  to  the  other  terms.  This  result  is  due  to 
the  boundary  layer  thickness  being  the  same  order  as  the  rotor  tip  radius 
and  much  larger  than  the  blade  chord.  These  gradients  are  not  as  large 
as  would  be  the  case  if  the  boundary  layer  thickness  would  be  much  le1.  * 
than  the  rotor  tip  radius.  However,  this  viscous  term  does  approach 
one  near  the  viscous  sublayer. 
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The  rate  of  change  of  the  streamwise  vortlcity  from  Equation  (10) 
simplifies  to 


3 rVl  2V  2VV  ^n'V 

3s'  ( V j WR'  w2  ~ w2 

Integrating  Equation  (11)  in  the  streamwise  direction  yields 


us'2  ws{  f Zu^.ds' 

4 

i 

1 

: fo)^,ds 

2 

t 2ft  ,u>,  ,ds 

I n'  b' 

W W = J WR' 

1 

T 

J 

w2 

L W 

J w2 

1 W 

(ID 


(12-a) 


where  'V*  W/R’ . 

Solving  Equation  (12-a)  for  ojs«  one  finds  that 


f 2un'  f 2fih,un'  ( 20  ,10.  , W 

“«2  " W2  J ds*  + W2  j ~2 ds*  - W2  j -—j ds’  + u)s,  ■ - 

1 1 W 1 W 1 


. (12-b) 


If  the  nonuniform  flow  through  the  rotor  was  such  that  the  relative 
streamlines  and  the  absolute  vortex  lines  lie  on  cylindrical  surfaces 

-X. 

(w^,*=0)  and  the  rotation  vector  (0)  is  parallel  to  the  axis  of  the 
cylindrical  surfaces,  then  Equation  (12-a)  would  become 


2io 


WR' 


ds' 


(13) 


This  relationship  is  commonly  used  in  most  cascade  flow  calculations. 
In  summary  for  a rotor  having  the  following  characteristics: 

1.  low  hub-to-tip  ratio, 

2.  rotor  operating  flow  coefficient  of  approximately  one 

3.  nonuniform  flow  over  the  span  of  the  blade. 
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the  following  two  equations  are  needed  to  estimate  the  passage  secondary 
flow: 


b) 

n 


W1  abj 

W ab, 


(7) 


2 

I 

1 

[3]  OTHER  APPROXIMATE  FORMS  FOR  THE  STREAMWISE  VORTICITY 

It  is  Important  in  secondary  flow  analysis  to  understand  the  flow 
field  in  order  to  determine  which  vorticity  terms  are  important.  Numerous 
analytical  and  semi-analytical  relationships  have  been  used  to  calculate 
the  streamwise  vorticity. 

Many  investigators  [12-16]  have  derived  expressions  for  cascade 
secondary  vorticity.  The  basic  analysis  is  due  to  Squire  and  Winter 
[12].  Their  expression  assumes  that  the  flow  is  incompressible  and 
inviscid.  Also,  they  assume  a constant  velocity  along  the  streamlines 
which  implies  that  the  normal  component  of  vorticity  remains  constant. 
Their  result  for  a rotating  frame  of  reference  is 


f 2w  . 

" W2  J wrT-  *■'  + W2 

1 


2W 


_W2J_^ds. 


ds' 


+ <*>s{ 


W„ 


(12- 


“s2  _ ws[  ° 2un,  e'  (14) 

where  e'  is  the  turning  angle  of  the  relative  streamline. 

Hawthorne  [13]  developed  a more  general  theory  for  the  secondary 
flow.  For  an  inviscid  and  incompressible  fluid  the  equation  written 
for  a rotating  frame  of  reference  would  be 
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This  equation  allows  only  for  a constant  axial  velocity  through  the  rotor 
and  a small  angle  of  turning. 

Marsh  [15]  has  extended  this  result  to  include  the  effect  of 
a change  in  axial  velocity  on  the  secondary  flow  at  exit  from  a cascade. 
It  is  shown  that  for  a row  of  inlet  guide  vanes  the  change  of  axial 
velocity  across  the  blade  row  has  a significant  effect  on  the  secondary 
vorticity.  His  expression  in  an  equivalent  rotating  frame  would  be 
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,x2 


2V 


X1 


sln201) 


x2 

- — tanB2~tan61 


tan62~tan8^ 


(17) 


Lakshminarayana  and  Horlock  [16]  have  extended  Equation  (15)  to  include 
the  effect  of  stream  surface  (Bernoulli  surface)  rotation  on  the  strength 
of  secondary  vorticity.  They  used  an  expression  by  Dean  [17]  in  which 
the  direction  and  magnitude  of  the  warping  of  real  stream  surfaces  is 
related  to  the  blade  geometry.  This  expression  in  terms  of  a rotating 
coordinate  system  is 


TT 

Y = 2 " 


2We|2  db' 


dW_e.2 


(18) 


where  the  streamlines  are  assumed  to  follow  the  curvature  of  the  blade 

camber  line.  As  can  be  seen  from  Equation  (18),  stream  surface  rotation 

becomes  important  when  the  incoming  velocity  gradients  (^T)  are  large 

db 

and  when  the  flow  is  turned  through  a large  angle  (e'). 

Using  Equation  (18)  in  Equation  (15)  gives  their  equation  for 
the  streamwise  vorticity  assuming  a constant  axial  velocity  and  the 
stagnation  pressure  varies  only  in  the  spanwise  direction.  Again,  the 
equivalent  expression  for  a rotating  reference  frame  is 


PL- 


U) 


2u)nJ 


W J cos^  cos82 


sin202  - sin261  &2~^1 

n a 


We'  “"i) 


®2~^l  82sln2^2-®lsin2^l 

t 


10 


f 
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83-g3 

83cos282-e3cos2B1  - - |(B2cos262-B1cos281) 


|(82sln282-B2sln2S1)  + |(sin2e2-sin281) 


(19) 


Dixon  [18]  developed  an  equation  which  can  be  applied  to  calculate 
the  secondary  flow  in  an  annular  cascade  having  a low  hub/tip  radius 
ratio.  In  his  analysis,  vortex  filaments  were  transported  by  a plane 
primary  flow  which  was  irrotational;  however,  the  stream  surfaces  of 
the  primary  flow  are  not  necessarily  at  the  same  radius  before  and  after 
the  cascade.  The  essential  result  for  the  outlet  streamwise  vorticity  is 


wS2 


^ COS^ 

— <DS'  ==■  + Idn. 

r2  S1  cos02  nl 


r^  sinf^ 
r2  sin82 


”l 

W„ 


tan£L  + 


r d 

J | _s 

2 T P2cosB2  J Wg 


(20) 


where  P2  is  the  distance  between  blades  at  exit. 


[ ds  „ r_ 

J W IT2 
s W„ 


(21) 


and 


F " n trl  Vrl  ~ r2  Vr2^ 


This  is  essentially  the  result  obtained  by  Hawthorne  and  Novak  [8]  for 
a plane  flow  in  a stationary  coordinate  system. 

Smith  [19]  has  presented  an  analysis  which  allows  secondary 
vorticity  to  be  calculated  for  mainstream  secondary  flows  using  the 
conventional  axisymmetric  solution  as  the  basic  flow,  i.e.,  the  flow 
is  assumed  to  remain  on  surfaces  of  revolution.  A major  assumption  of 


(22) 
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Smith's  analysis  Is  that  the  distortion  of  those  stream  surfaces  Is  not 
large.  A relation  for  the  secondary  streamwise  vorticity  In  the  flow  at 
the  trailing  edge  plane  was  achieved  with  aid  of  the  vortex  laws.  The 
expression  obtained  is 


where  a is  the  distance  between  the  exit  streamlines,  T,,.  is  the  blade 

VA 

circulation  in  the  actual  flow,  i.e.,  primary  and  secondary  flow,  and 
Ty  is  the  blade  circulation  in  the  primary  flow. 

It  is  important  to  note  that  there  are  differences  in  the  definition 
of  the  primary  flow  between  Smith's  theory  and  other  theories.  However, 
Horlock  [5]  shows  that  all  approaches  lead  to  basically  the  same  cascade 
secondary  flow  theory. 


[4]  SUMMARY 

Equations  for  the  development  of  streamwise  secondary  vorticity  in  a 
blade  passage  for  rotating  systems  are  derived  from  vector  equations  for 
vorticity  using  intrinsic  coordinates.  Comparisons  can  be  made  between  the 
secondary  vorticity  equations  for  a rotor  operating  at  a flow  coefficient 
near  one  and  those  which  are  commonly  used  in  cascades. 

The  equations  for  a rotor  operating  at  a flow  coefficient  near  one 
show  a coupling  of  equations  for  the  vorticity  components  parallel  and 
normal  to  the  streamlines.  In  order  to  solve  this  set,  the  equations 
were  nondimensionalized  with  respect  to  the  characteristic  parameters 
for  this  three-dimensional  flow.  The  resulting  equation  for  the  stream- 
wise  vorticity  allows  for  convergence  of  the  streamlines  in  the  radial 


J 
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direction,  for  the  effects  of  rotation,  and  for  the  effects  of  streamline 
curvature. 

Streamwise  component  of  vorticity  Is  produced  In  a blade  passage 
through  deflection  of  the  relative  flow  having  a normal  vorticity 
component.  Additional  effects  on  the  streamwise  vorticity  occur  due  to 
rotor  rotation  and  streamline  curvature.  It  is  interesting  to  note  that 
both  a relative  velocity  increase  through  the  rotor  and  stream  surface 
rotation  decreases  the  amount  of  vorticity  generated  by  the  turning  of  the 
normal  vorticity  vector  through  the  rotor. 
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Figure  3 - Notation  for  Rotating  Coordinate  System 
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BERNOULLI 


Figure  4 - Flow  Through  a Blade  Passage 


